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Abstract 

The drying of food can extend the shelf life of food, reduce transportation and storage costs. Fick's second 
law is commonly used to evaluate the mass data in the drying process in a standard way and is based on many 
assumptions. Understanding the meaning of mass transfer in products can improve the drying process and 
product quality. Computational fluid dynamics (CFD) models fluid flow situations utilizing powerful computer 
and applied mathematics in order to predict mass transfer in industrial processes. The aim of this research was 
numerical study of the drying behavior of pistachio nut using CFD method and evaluating the numerical results 
in the bed condition of fluid, semi fluid and fix bed as well as air temperatures of 90, 75, 60 and 45°C. During 
drying using computational fluid dynamic and the Fluent CFD code, the external flow and temperature fields 
around the cylindrical object (7.5x 10 millimeter) will be predicted in the numerical analysis. A laboratory 
fluid bed dryer was used for drying experiments. The main parts of the dryer are forward radial fan, drying 
chamber, electrical heater, inverter, temperature controller. The dryer attachment tools are input and output 
temperature sensors, anemometer and computer. The numerical part was verified and juxtaposed with the 
experimental data. The numerical solution result at 60, 75 and 90°C were so close to experimental results 
except for air temperature of 45°C. Mean absolute error in fix bed, at 60, 75 and 90°C were 0.2123, 0.1257 
and 0.0337 which were lower than 45°C temperature and R? values for these temperatures were 0.9903, 0.9705 
and 0.9807, respectively. As the temperature decreased, the values of Ew, and X’ increased in all bed 
conditions. The average value of R? for all applied bed conditions was 0.9850. This value showed high 
correlation between experimental and numerical results. 


Keywords: Semi fluid bed drying, Fluid bed drying, Fix bed drying, Experimental results, Modeling, 
Numerical solution. 


Introduction 

Nuts are rich in nutrients, fiber, protein, 
phytosterols and antioxidants (Fantino et al., 
2020). The edible seeds harvested from the 
pistachio tree are pistachios. Among the 
different species of the genus Pistachio, 
pistachio (Pistacia Vera) is the only known 
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commercially edible nut. Its main cultivation 
areas include the Mediterranean region. It is not 
affected by drought and salinity (Noguera 
Artiaga et al., 2020). For consumers to consider 
it an edible nut, the fruit of the Pistacia Vera 
variety must be large enough (Shokraii and 
Esen, 1998). Pistachios are used as ingredients 
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in the confectionery, snacks, ice cream and 
pastry industries. Pistachios were originally 
grown in the United States, the Middle East, 
and especially in Iran (Kouchakzadeh and 
Shafeei, 2010).The quality of pistachios 
depends on proper harvesting and post-harvest 
processing. Drying is very important in the 
marketing of the final product. The optimal 
moisture content for storing pistachios is 7%- 
9% (w.b.) (USDA Std., 1990). 

Studies have shown that drying methods and 
different changes are key factors in determining 
the quality of dried food products (Naidu et al., 
2016). In other words, the deterioration and 
damage of certain characteristics, including 
their color, structure, aromatic compounds and 
nutrients, may lead to a decline in product 
quality (Izli and Isik, 2015). The most common 
method used to dry pistachios is the hot air 
drying method. Fluidized bed drying is a 
profitable drying method for pistachios. When 
passing through the pistachio bed, a fixed bed 
is formed when a low-velocity airflow flows 
upward. When the highest air velocity is 
applied, the whole pistachio starts to float, 
which is called minimum _ fluidization 
(AmiriChayjan et al., 2012). The semi- 
fluidized bed condition occurs at the maximum 
static pressure drop. When the inflow airflow 
increases, a bubbling fluidized bed is created, 
causing the mixing of pistachios (Kuni and 
Levenspiel, 1991; Kaveh and AmiriChayjan, 
2015). 

In the drying mechanism of a system, the 
phenomenon of simultaneous transmission 
(energy, momentum and mass) will occur. 
Drying can extend the shelf life of food, reduce 
transportation and storage costs, and develop 
new consumption methods (Lindsay Rojas and 
Augusto, 2018). Fick's second law (Fick, 1855) 
is commonly used to evaluate the mass data 
process in the drying process in a standard way 
and is based on many assumptions (Lindsay 
Rojas and Augusto, 2018). Understanding the 
meaning of mass transfer in products can 
improve the drying process and product quality. 
External factors include temperature, relative 
air humidity, and air velocity. However, 
internal factors include factors such as the 


density and permeability of the material, as well 
as thermophysical properties (Kaya ef al., 
2006). 

Compared with experimental research, 
numerical simulation is a common method 
applied for drying process analysis due to lower 
cost and time (Kaya et al., 2006). Numerical 
methods for describing food drying 
mechanisms can provide useful information to 
help understand temperature and humidity 
requirements more clearly (Haghighi ef al., 
1990; Rafiee and Kashaninejad, 2005; Rafiee et 
al., 2005). In computational fluid dynamics 
(CFD), numerical methods are often used to 
approximate the equation that governs fluid 
dynamics in the target fluid region. It uses 
numerical methods to predict chemical 
reactions, mass transfer, heat transfer, fluid 
flow, and related phenomena by solving 
mathematical equations that control these 
processes (Puma Chandra, 2017). CFD 
technology can also help solve complex 
transmission phenomena and make drying 
process more cost and time-effective (Kaya et 
al., 2006; Kaya et al., 2008a). 

Many researchers studied numerical 
modeling of mass transfer for a wide variety of 
agricultural and food products. These include 
hazelnut (Topuz ef al., 2004), garlic 
(AbbasiSouraki and Mowla, 2008), kiwi fruit 
(kaya et al., 2008b), apple slices (Mabrouk et 
al., 2012), soybean meal (Silva et al., 2012), 
shrinkable products such as lentils (Carmo and 
Lima, 2004), fruits and vegetables (Kowalski 
and Mierzwa, 2013). However, no studies have 
specifically examined numerical modeling for 
pistachio during convective drying of fluid, 
semi fluid and fix bed. 

Researchers around the world have studied 
the drying of various products using CFD 
simulation (Malekjani and Jafari, 2018; 
Demissie et al., 2019).Studies have also been 
performed on the simulation of pistachio nuts 
drying by finite element (Rafiee et al., 2007; 
Rafiee et al., 2009). No study has been done on 
simulating the drying of pistachio nuts using 
CDF simulation. The aim of this research was 
numerical study of the drying behavior of 
pistachio nut using CFD method and e 
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evaluating the numerical results in the bed 
condition of fluid, semi fluid and fix bed as well 
as air temperatures of 90, 75, 60 and 45°C. 
During drying process using computational 
fluid dynamic and the Fluent CFD code, the 
external flow and temperature fields around the 
cylindrical object (7.5x 10 millimeter) were 
predicted in the numerical analysis. The 
numerical part was verified and_ the 
experimental data relevant to the method of 
thin-layer drying were juxtaposed. 


Materials and method 

Experimental setup and procedure 

Ohadi cultivar of pistachio was selected for 
conducting the study. The initial moisture 
content of the samples collected was 50.3+0.2% 
(d.b.), at 130°C and 24 h using the oven process 
(AOAC, 1995). 

To indicate bed conditions (air velocities) in 
tests, the increasing air velocity against 
pressure drop were recorded and then plotted as 
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Fig. 1. An estimating and recording unit for 
differential digital manometer and a vane type 
advanced anemometer (StandardST-8897) was 
utilized to get the fluidization curve. To 
accomplish the net air pressure drop across the 
pistachio bed, from the start, total static 
pressure drop because of pistachio column and 
bed plate was estimated, and afterward was 
deducted from air pressure drop because of 
empty chamber. 

The greatest value of static pressure drop 
against a specific air velocity in fluidization 
curve is characterized as least fluidization point 
or semi fluidized bed (Kunii and Levenspiel, 
1991). Fluidization tests were acted in three 
recreates for thin layer drying of pistachio tests 
with around 100 g load. In the wake of getting 
the semi fluid bed point (air velocity about 2.6 
ms"'), two focuses with air velocities of 1.6 
ms! and 4.1 ms‘! were chosen as fix and fluid 
bed conditions, respectively. 


28 3 38 4 45 =§& 
Airvelocity (m/s) 
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Fig. 1. Pistachio nuts' fluidization curve 


The drying process was ended, when the 
moisture content of the samples reached an 
average of about 9 percent (d.b.) (Amirichayjan 
et al., 2012). The air relative humidity and air 
temperature were 21%-33% and 28-32°C, 
respectively. A laboratory fluid bed dryer was 
used for drying experiments. This apparatus 
was fabricated in agricultural machinery 
engineering of Bu-Ali Sina University (Figure 
2). The main parts of the dryer are forward 


radial fan, drying chamber, electrical heater, 
inverter, and temperature controller. The dryer 
attachments tools are input and _ output 
temperature sensors, anemometer and 
computer. 45, 60, 75 and 90°C air temperature 
levels and fix, semi fluid and fluid bed 
conditions were used as input variables in the 
experiments. Each experiment was replicated 
three times (Amirichayjan et al., 2012). 
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In each drying test, around 100 g pistachio 
sample was stacked in drying chamber. 
Pistachio nuts were extended in a thin layer 
form in the drying chamber and the test was 
begun. The ambient air temperature, input and 
output air temperature, air velocity, air relative 
humidity and sample weight were consistently 
checked and recorded during drying process. A 


Fan and electro motor 


Hygrometer / Etectrical heater 


Chassis Electrical panel 


digital balance (AND GF-6000) with 0.01 g 
accuracy was used to online weighing the nuts 
during the drying experiments. Moisture 
content of the pistachio samples in each drying 
run was calculated based on the initial and final 
moisture content of the samples and initial 
samples mass. 


Drying chamber Air velocity recorder 


Outlet air temperature recorder 
iF 


Input air temperature recorder 


Fig. 2. Laboratory Scale Schematic Diagram of the Fluidized Bed Dryer (Golpour et al, 2021) 


Modeling 

Partial differential equations that govern the 
forced convection motion of a fluid bed drying 
in a two-dimensional geometry are included in 
the energy, mass and momentum conservation 
equations. In some cases that are not complex, 
physical and thermal properties are considered 
constant. For a two-dimensional cylindrical 
problem, considering the flow incompressible, 
the Navier-Stokes equations are described in 
their most general form: 


The equation of mass conservation (i.e. 
continuity) is (Norton and Sun, 2006): 
(1) 


Ee Gea a =0 
r Or Oz 


where the radial coordinate is , and the axial 
coordinate is z. 

The momentum equations are (Norton and 
Sun, 2007): 
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where p represents density (kg/m), yw 
denotes the dynamic viscosity (Pa.s), pis the 


Uu, | (2) 


(3) 


pressure (Pa) and ushows the velocity in x- 
direction (m/s). 
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The energy equation is (Ferziger and Peric, 
20002): 


10, OT 0 OT 
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which 7 shows the air temperature (K), fis 
the drying time (h), vrepresents the velocity in 
y-direction (m/s) and a@ denotes the thermal 
diffusivity (m7/s). 

The Fluent V6.3.26 CFD package based on 
the finite volume method has been used for 
transforming and solving these equations. The 
boundary conditions assumed are as follows: 

1) For velocity, no-slip conditions 

2) The constant surface temperature of 
the drying material. 

The following governing equations related 
to two-dimensional heat and moisture transfer 
can be written as follows, taking into account 
the above assumptions (Anderson, 1992): 
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where C,, represents the constant pressure 


specific heat (J/kg K), k stands for the thermal 
conductivity (W/m K), Dis the moisture 
diffusivity (m?/s) and W denotes the moisture 
content (kg/kg, d.b.). 


In both instances, the following initial and 
boundary conditions were used (Anderson, 
1992): 


T(r, z,0) =T, and M(r,z,0) =M, 
Where M,is the initial moisture content 
(kg/kg, d.b.), nis normal to surface, s 


represents the surface coordinate, h denotes the 
heat transfer coefficient (W/m? K) and h, shows 


the mass transfer coefficient (m/s). 
The spectrum of coefficients associated with 
convective heat and mass transfer (handh,,) on 


the surface of the material was assumed. For 


! 
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At 
Or: 

(5) 


solving the mass transfer equations under the 
initial and boundary conditions of interest, the 
finite difference method was applied. 


k Aa por=A(T, a T,,) 
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Oz 20M (M, M.,) 
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Various researchers have used the following 
correlations for calculating parameters for the 
pistachio, as follows: 

For bulk density (Hsu ef al., 
p= 4394+5.003M (R? =0.959 ) 


1991); 


For thermal diffusivity (Hsu er al., 1991); 
a@=51.1x10° —0.568x 10M (R? =0.983) 


For thermal conductivity (Hsu et al., 1991); 
k =0.0866+0.2817x10°M (R? =0.963) 

For constant-pressure specific heat (Hsu et 
al., 1991); C, =1074+27.79M (R* =0.920) 


For moisture diffusivity (AmiriChayjan et 
al.,2012); Dy =4x10°m"/s 
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Finally, Four indices such as correlation 
coefficient (R*), absolute error (£,,,,), mean 


squared error (MSE) and chi-square (y”) were 
used as the goodness of fit and agreement 
between experimental results and numerical 
solutions. These indices are as_ follow 
(AmiriChayjan et al., 2012): 


N (9) 
> [Meco: a M ans 
R =] i=l 
N = M umi 
2 M ons = NO 


M 


numi 


1 “ exp.i a 
Pe = 2 » M 


expJ 


| (10) 


which M represents the experimental 


moisture ratio of idata, M 


exp.i 


denotes the 


numi 


numerical moisture ratio of i” data and N 
shows the number of observations. 


1x , 
MSE = 5) Ui-y) 
t=1 


where N_ stands for the number of 
observations, f; represents the numerical value 
of i data and y; shows the experimental value 


of i” data. 
N 


Ais »: (0; — E;)? 
i E; 
l=1 


In which N denotes the number of 


observations, O; is the experimental value of 7 a 


data and E; is the numerical value of i” data. 


Results and discussion 
Comparison of experimental and numerical 
results 


Drying curves of pistachio nuts in all the 
three bed conditions and _ temperature 
conditions using both experimental and 
numerical methods are presented in Figs. 3-5. 

After drying the pistachio nuts to nearly 9% 
(d.b.), the simulated moisture contents at 45°C 
were found to be higher compared to the 
measured values (Fig. 3A). Results related to 
the simulated moisture curve and to the 
experimental variation of moisture at an air 
temperature of 60°C are shown (Figure 3B). 
Based on this result, the simulated values within 
the range of 6000 and 10000 seconds were 
slightly higher than the measured values. 
Because the properties of the product are 
expressed as averages, so for half of the 
conditions, the forecast is done properly 
(Makarichian et al., 2021; Rashidi et al., 2021). 

At 75°C the simulated values between 8500 
and 30000 seconds were slightly lower and 
those between 47000 and 66000 seconds were 
higher compared to the measured values (Fig. 
3C) and at 90°C the values simulated were 
lower than the values measured (Fig. 3D). 

Fig. 4A till Fig. 4D show the simulated 
moisture contents for semi fluid bed at 45°C 
and 90°C, respectively. They had the same 
trend in contrast with fix bed. The reason for 
this condition can be the average amount of 
physical and thermal properties 

At 60°C the simulated values between 60000 
and 10800 seconds were slightly higher than the 
measured values (Fig. 4B). Fig. 4C showed that 
the simulated values from 48000 to 66000 
seconds were slightly higher than the measured 
values (at 75°C). 

The same trends were found for fluid, fix and 
semi fluid bed at 45°C and 90°C (Fig. 5A and 
Fig. 5D), also the curve trends for fluid bed was 
the same against semi fluid bed at 60°C and 
75°C (Fig. 5B and Fig. 5C). 

According to Figs. 3- 5, following an 
increase in air temperature, drying time 
decreased. 
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Fig. 3: Measured and simulated moisture contents at fix bed (velocity=1.6 m/s). 
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Fig. 4. Measured and simulated moisture contents at semi fluid bed (velocity=2.6 m/s). 
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Fig. 5. Measured and simulated moisture contents at fluid bed (velocity=4.1 m/s). 
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The rate of energy exerted to the pistachio 
sample bed and then the drying rate increased 
following increased the temperature being 
applied to the pistachio bed. Increased 
temperature allows the rate of transfer of heat to 
the product to be increased. Air velocity also 
has no major effect on the rate of drying. An 
increase in the coefficients of convective heat 
and mass transfer between the drying air and the 
pistachio nuts caused an increase in the rate of 
drying and decreased the drying time. 


Similar results were reported by other 
studies on drying other products including 
hazelnuts (Topuz et al., 2004), eggplant 
(Akpinar and Bicer, 2005), peach (Kingsly et 
al., 2007), plum (Goyal et al., 2007), berberis 
fruit (Aghbashlo et al., 2007), corn (Vukié et 
al., 2015), ginger (Parlak, 2015), chilean berry 
(Quispe-Fuentes ef al., 2016) and turnip (Kaveh 
and AmiriChayjan, 2016). 

Table 1 shows the values related to R’, E.,,,, 


MSE and y2, which were calculated for all 
examined bed conditions. The highest and the 


lowest values of rR? and E,,,, confirm the 


highest agreement between 
results and numerical solutions. 
Figure 6 A-C shows the comparison between 
the numerical solution and experimental results 
on each bed condition. 
According to Fig 6, it seems that the 
numerical solution result at 60, 75 and 90°C 


experimental 


were so close to experimental results except for 
air temperature of 45°C. The values of each 
parameter in Table 1 proved this result. In fix 
bed, the value of E,,,at 60, 75 and 90°C were 


0.2123, 0.1257 and 0.0337 that were lower than 
the air temperature of 45°C (0.2595) and R? 
values for these temperatures were 0.9903, 
0.9705 and 0.9807. These trends were repeated 
for the bed conditions of semi fluid and fluid 
bed with different values. 

The average values of R* for all applied bed 
conditions was calculated 0.9850. This value 
showed the high agreement between 
experimental results and numerical solutions. 

According to the values of statistical indices 
for different conditions, the errors created by 
using the model type and fixed numbers of 
physical and thermal properties are not 
significant and the results of this study can be 
used to predict the kinetics of dried pistachios 
at any temperature and used air velocity in the 
study area. The advantage of this numerical 
method is its accuracy and comprehensiveness. 
However, to predict the drying kinetics of 
pistachios with the help of experimental and 
semi experimental models, it will be necessary 
to adjust the model coefficients (Azharul Karim 
and Hawlader, 2005; Shahbazi and Rahmati, 
2013; Castro et al., 2018; Nguyen ef al., 2021). 


Table 1- Values of R?, Eass, MSE and X’ for all applied bed conditions 


Bed condition Temperature (°C) R? Eabs MSE xX? 
Fix bed 0.9955 0.2595 0.0017. 03 
(1.6 m/s) 60 0.9903 0.2123 0.0010 0.13 

75 0.9705 0.1257 0.0005 0.059 
90 0.9807 0.0337 0.0004 0.018 
Semi fluid bed 45 0.9944 0.1850 0.0008 0.15 
(2.6 m/s) 60 0.9879 0.1706 0.0007 0.098 
75 0.9715 0.1475 0.0005 0.053 
90 0.9856 0.0120 0.0003 0.017 
Fluid bed 45 0.9942 0.2908 0.0022 0.34 
(4.1 m/s) 60 0.9911 0.1746 0.0008 0.098 
75 0.9808 0.1278 0.0004 0.045 
90 0.9782 0.0662 0.0007 0.029 


114 Iranian Food Science and Technology Research journal, Vol. 17, No. 6, Feb. Mrch. 2022 


- Numerical and Experimental solution at fix bed 
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Fig. 6. Experimental results and the numerical solution of drying kinetic of pistachio nuts in bed conditions 
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Conclusions 

The results of this research work are as 
follows: 

1) Comparison between the experimental 
data with the results obtained in the numerical 
solution with Fluent CFD code showed that 
there is a high agreement between them. 

2) — Infix bed (1.6 ms”), the lowest values 
of E,,,, (0.0337) and 72 (0.0004) were related 


to the temperature of 90°C. The R? value of this 
temperature was 0.9807, which shows the 
highest agreement between experimental data 
and numerical solution compared to other 
examined temperatures. 

3) The lowest values of v2 and E_,, in 


the semi fluid bed (2.6 ms!) were 0.017 and 
0.0120 and belong to the temperature of 90°C 
with R’= 0.9856. In semi fluid bed, the highest 


data and numerical solution for 90°C. After 
90°C, temperatures of 75, 60 and 45°C showed 
the greatest agreement between experimental 
and numerical results, respectively. 

4) The temperature of 90°C showed the 
highest agreement in the fluid bed (4.1 m s-1). 
The values of R*, v2 and E,,. for this 


temperature were 0.9782, 0.029 and 0.0662, 
respectively. 

5) The average value of R? for all applied 
bed conditions was calculated 0.9850. This 
value showed the high agreement between 
experimental results and numerical solutions. 

6) The highest agreement in all bed 
conditions was related to temperature of 90°C. 

7) As the temperature decreased, the 
values of E.,, and v2 increased in all bed 


AY 


conditions. 


agreement was obtained between experimental 


References 

AbbasiSouraki B., Mowla D., 2008, Experimental and Theoretical Investigation of Drying Behavior 
of Garlic in an Inert Medium Fluidized Bed Assisted by Microwave. Journal of Food Engineering, 
88, 438— 449. 

Aghbashlo M., Kianmehr M.H., Samimi-Akhijahani H., 2008, Influence of Drying Conditions on the 
Effective Moisture Diffusivity, Energy of Activation and Energy Consumption During the Thin- 
Layer Drying of Beriberi Fruit (Berberidaceae). Energy Conversion and Management, 49, 2865— 
2871. 

Akpinar E.K., Bicer Y., 2005, Modeling of the Drying of Eggplants in Thin Layer. International 
Journal of Food Science and Technology, 40, 273— 281. 

AmiriChayjan R., Alizade H.H.A., Shadidi B., 2012, Modeling of Some Pistachio Drying 
Characteristics in Fix, Semi Fluid and Fluid Bed Dryer. Agricultural Engineering International: 
CIGR Journal, 14, 143- 154. 

Anderson J.D., 1992, Governing Equations of Fluid Dynamics. In: Wendt J.F. (eds) 
Computational Fluid Dynamics. Springer, Berlin, Heidelberg. 

AOAC, 1995, Official Methods of Analysis, Association of Official Analytical Chemists, 
Washington, DC, USA. 

Azharul Karim Md., Hawlader M.N.A., 2005, Mathematical modelling and experimental 
investigation of tropical fruits drying, International Journal of Heat and Mass Transfer, 48, 4914- 
4925. 

Carmo J.E.F., Lima A.G.B.D., 2004, Drying of Lentils Including Shrinkage: A Numerical 
Simulation. Drying. Proceedings of the 14th International Drying Symposium, Sao Paulo, Brazil, 
22-25, 510-517. 

Castro A.M., Mayorga E.Y., Moreno F.L., 2018, Mathematical modelling of convective drying of 
fruits: A review, Journal of Food Engineering, 223, 152- 167, 

Demissie P., Hayelom M., Kassaye A., Hailesilassie A., Gebrehiwot M., Vanierschot M., 2019, 
Design, development and CFD modeling of indirect solar food dryer. Energy Procedia, 158, 1128- 
1134. 


116 Iranian Food Science and Technology Research journal, Vol. 17, No. 6, Feb. Mrch. 2022 


Fantino M., Bichard C., Mistretta F., Bellisle F., 2020, Daily consumption of pistachios over 12 weeks 
improves dietary profile without increasing body weight in healthy women: A randomized 
controlled intervention. Appetite, 144, 104483. 

Ferziger J.F., Peric M., 2002, Computational methods for fluid dynamics (3rd ed.). New York: 
Springers publications. 

Golpour I., Guine R.P.F., Poncet S., Golpour H., Amirichayjan R., AmiriParian J., 2021, Evaluating 
the heat and mass transfer effective coefficient during the conductive drying process of paddy. 
Journal of Food Process Engineering, 13771. 

Goyal R.K., Kingsly A.R.P., Manikantan M.R., Ilyas S.M., 2007, Mathematical Modeling of Thin 
Layer Drying Kinetics of Plum in a Tunnel Dryer. Journal of Food Engineering, 79, 176— 180. 
Haghighi K., Irudayaraj J., Stroshine R.L., Sokhansanj S., 1990, Grain Kernel Drying Simulation 

Using the Finite Element Method. Transactions of the ASAE, 33, 1957— 1965. 

Hsu M.H., Mannapruma D., Singh R.P., 1991, Physical and Thermal Properties of Pistachios. Journal 
of Agricultural Engineering Research, 49, 311- 321. 

Izli N., Isik E., 2015, Color and Microstructure Properties of Tomatoes Dried by Microwave, 
Convective, and Microwave-Convective Methods. International Journal of Food Properties, 18, 
241- 249. 

Kashaninejad M., Tabil L.G., Mortazavi A., Safekordi A., 2003, Effect of Drying Methods on Quality 
of Pistachio Nuts. Dry Technology, 21, 821— 838. 

Kaveh M., AmiriChayjan R., 2015, Mathematical and Neural Network Modeling of Terebinth Fruit 
Under Fluidized Bed Drying. Research in Agricultural Engineering, 61(2), 55— 65. 

Kaveh M., AmiriChayjan R., 2016, Modeling Thin-Layer Drying of Turnip Slices under Semi- 
Industrial Continuous Band Dryer. Journal of Food Processing and Preservation, 41(2), 1- 14. 
Kaya A., Aydin O., Dincer I., 2006. Numerical Modeling of Heat and Mass Transfer during Forced 
Convection Drying of Rectangular Moist Objects. International Journal of Heat and Mass 

Transfer, 49, 3094— 3103. 

Kaya A., Aydin O., Dincer I., 2008a, Heat and Mass Transfer Modeling of Recirculating Flows 
During Air Drying of Moist Objects for Various Dryer Configurations. Numerical Heat Transfer 
Part A-Applications, 53: 18— 34. 

Kaya A., Aydin O., Dincer I., 2008b, Experimental and Numerical Investigation of Heat and Mass 
Transfer During Drying of Hayward Kiwi Fruits (ActinidiaDeliciosa Planch). Journal of Food 
Engineering, 88, 323-330. 

Kingsly A.R.P., Goyal R.K., Manikantan M.R., Ilyas S.M., 2007, Effects of Pretreatments and 
Drying Air Temperature on Drying Behavior of Peach Slice. International Journal of Food 
Science and Technology, 42: 65- 69. 

Kunii D., Levenspiel O., 1991, Fluidisation Engineering, Butterworth-Heinemann. 

Kouchakzadeh A., Shafeei S., 2010, Modeling of Microwave-Convective Drying of Pistachios. 
Energy Conversion and Management, 51: 2012-2015. 

Kowalski S.J., Mierzwa D., 2013, Numerical Dnalysis of Drying Kinetics for Shrinkable Products 
Such as Fruits and Vegetables. Journal of Food Engineering, 114, 52- 59. 

Lindsay Rojas M., Augusto P., 2018, Microstructure elements affect the mass transfer in foods: The 
case of convective drying and rehydration of pumpkin. LWT - Food Science and Technology 93, 
102- 108. 

Mabrouk S.B., Benali E., Oueslati H., 2012, Experimental Study and Numerical Modelling of Drying 
Characteristics of Apple Slices. Food and Bio products Processing, 90, 719-728. 

Makarichian A., Amiri Chayjan R., Ahmadi E., Mohtasebi S.S., 2021, Assessment the influence of 
different drying methods and pre-storage periods on garlic (Allium Sativum L.) aroma using 
electronic nose. Food and Bioproducts Processing, 127, 198- 211. 


Shadidi and Amiri Chayjan/ Mass transfer simulation of pistachio nuts using .... 117 


Malekjani N., Jafari S.M., 2018, Simulation of food drying processes by Computational Fluid 
Dynamics (CFD); recent advances and approaches. Trends in Food Science & Technology, 78, 
206- 223. 

Naidu M.M., Vedhashree M., Satapathy P., Khanum H., Ramsamy R., Hebbar H.U., 2016, Effect of 
Drying Methods on the Quality Characteristics of Dill (Anethumgra-Veolens) Greens. Food 
Chemistry, 192, 849- 856. 

Nguyen T.-D., Nguyen-Quang T., Venkatadri U., Diallo C., Adams M., 2021, Mathematical 
Programming Models for Fresh Fruit Supply Chain Optimization: A Review of the Literature and 
Emerging Trends. AgriEngineering, 3, 519-541. 

Noguera-Artiaga L., Sdnchez-Bravo P., Pérez-Lépez D., Szumny A., Calin-Sanchez A., Burgos- 
Hernandez A., Carbonell-Barrachina A., 2020, Volatile, Sensory and Functional Properties of 
HydroSOS Pistachios. Foods, 92(2), 158. 

Norton T., Sun D.W., 2006, Computational fluid dynamics (CFD) an effective design and analysis 
tool for the food industry: a review. Trends in Food Science and Technology, 17(11), 600— 620. 

Norton T., Sun D.W, 2007, An overview of CFD applications in the food industry, chapter 1. In D.- 
W. Sun (Ed.), Computational fluid dynamics in food processing (pp. 1-43). Boca Raton: CRC 
Press. 

Parlak N., 2015, Fluidized Bed Drying Characteristics and Modeling of Ginger (ZingiberOfficinale) 
Slices. Heat and Mass Transfer, 51 (8), 1085-1095. 

Puma Chandra B., 2017, Introduction to computational fluid dynamics. International Journal of 
Science and Computing, 3(2), 177- 120. 

Quispe-Fuentes I., Vega-Galvez A., Vasquez V., Uribe E., Astudillo S., 2016, Mathematical 
Modeling and Quality Properties of a Dehydrated Native Chilean Berry. Journal of Food Process 
Engineering, 32, 56- 69. 

Rafiee Sh., Kashaninejad Mahdi., Keyhani A., Jafari A, 2009, Pistachio Nut (Ohadi Variety) Mass 
Transfer Simulation during Process of Drying Using Finite Element Method. Journal of 
Agriculture and Science Technology, 11, 137- 146. 

Rafiee S., Kashaninejad M., 2005.Transient Moisture Gradients in Pistachio Nut with Finite Element 
Model during High Temperature Drying. /V International Symposium on Pistachio and Almonds, 
Book of Abstracts,Tehran, Iran. 22- 25. 

Rafiee S., Jafari, A., Kashaninejad M., Omid M., 2007, Experimental and numerical investigations 
of moisture diffusion in pistachio nuts during drying with high temperature and low relative 
humidity. International Journal of Agriculture and Biology, 9, 412- 415. 

Rafiee S., Kashaninejad M., Tabatabaeefar A., 2005, Transient Moisture Gradients in Wheat (Tagan) 
Kernel with Finite Element Model. Asia Pacific Drying Conference. December 13-15, Kolkata, 
India, 732- 740. 

Rashidi M., Amiri Chayjan R., Ghasemi A., Ershadi A., 2021, Tomato tablet drying enhancement by 
intervention of infrared - A response surface strategy for experimental design and optimization. 
Biosystems Engineering, 208, 199- 212. 

Shokraii E.H., Esen A., 1998, Composition, Solubility and Electrophoretic Patterns of Protein 
Isolated from Kerman Pistachio Nuts (Pistaciavera L.). Journal of Agricultural and Food 
Chemistry, 36, 425-429. 

Shahbazi F., Rahmati S., 2013, Mass modeling of fig (Ficus carica L.) fruit with some physical 
characteristics. Food Sci Nutr. 1(2), 125- 129. 

Silva F.R.G.D., Souza M.D., Costa A.M.D.S.D., Jorge L.M.D.M., Paraiso P.B., 2012, Experimental 
and Numerical Analysis of Soybean Meal Drying in Fluidized Bed. Powder Technology, 229, 61— 
70. 

Topuz A., Gur M., ZaferGul M., 2004, An Experimental and Numerical Study of Fluidized Bed 
Drying of Hazelnuts. Applied Thermal Engineering, 24, 1535— 1547. 


118 Iranian Food Science and Technology Research journal, Vol. 17, No. 6, Feb. Mrch. 2022 


USDA Standard. 1990. US Standard for grades of pistachio nuts, Agri Marketing Service, 
Washington DC, USA. 

Vukié M., Janevski J., Vuckovié G., StojanovicB., Petrovié A., 2015, Experimental Investigation of 
the Drying Kinetics of Corn in a Packed and Fluidized Bed. Iranian Journal of Chemistry and 
Chemical Engineering, 34(3), 43-49. 


Nomenclatures 
a thermal diffusivity (m7/s) P pressure (Pa) 
C, constant pressure specific heat (J/kg K) r radial coordinate 
D moisture diffusivity (m/s) KY surface coordinate 
Be absolute error T air temperature (K) 
h heat transfer coefficient (W/m? K) t Drying time 
h,, mass transfer coefficient (m/s) u velocity in x direction (m/s) 
k thermal conductivity (W/m K) v velocity in y direction (m/s) 
M xpi experimental moisture ratio W moisture content (kg/kg, db) 
M vu numerical moisture ratio Zz axial coordinate 
M, initial moisture content (kg/kg, db) pP density (kg/m?) 
N number of observations i dynamic viscosity (Pa.s) 
n normal to surface 


correlation coefficient 


a 
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